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Search for CP violation in the decay — > D*^D^ 
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Abstract 

We report a search for CP-violating asymmetry in D*^D^ decays. The analysis employs 

two methods of reconstruction: full and partial. In the full reconstruction method all daughter 
particles of the B^ are required to be detected; the partial reconstruction technique requires a 
fully reconstructed and only a slow pion from the D*'^ -^"^siow decay. From a fit to the 
distribution of the time interval corresponding to the distance between two B meson decay points 
we calculate the CP- violating parameters and find the significance of nonzero CP asymmetry to 
be 2.7 standard deviations. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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In the Standard Model (SM), CP violation arises from the Kobayashi-Maskawa (KM) 
phase ^ in the weak interaction quark-mixing matrix. Comparisons between SM expec- 
tations and measurements in various modes are important to test the KM model. The 
D*'^D^ modes are of particular interest since large CP violation effects are expected 
in these decays j^l- Although the D*^D^ final states are not CP eigenstates, they can be 
produced in the decays of both and B^ with comparable amplitudes; the interference 
between amplitudes of the direct transition and that via BB mixing results in CP violation. 
These decays are dominated by the tree h — ccd transition, thus CP violation measurements 
are sensitive to the angle 0i. However, the h ^ d penguin diagram also contributes to this 
final state and contains a different weak phase. Therefore this contribution results in both 
direct CP violation and a deviation of the mixing-induced CP violation parameter from 
sin 201. The Cabibbo suppressed decays B^ D*^D^ were first observed by Belle jsj, and 
have been confirmed by BABAR 0]. 

The probability for a B meson to decay to D*^D^ at time At can be expressed in terms 
of five parameters. A, S± and C±: 

~\At\/TgO 

P±,^(At) = (1 ± A)— {1 + q[S± sin(Am,,At) - C± cos(AmrfAt)]}. (1) 

OTbO 

Here the +(— ) sign represents the D*~^D~ {D*~D^) final state, and the 6-fiavor charge 
q = +1(— 1) when the tagging B meson is a B^{B^). The time-integrated asymmetry A 
between the rates to D*^D~ and D*~ is defined as 

Nd*+d- + Nf)*-jj+ 

In the case of negligible penguin contributions , the parameters S± can be related to 
the weak phase difference (sin20i in the SM), the strong phase difference (5) and the ratio 
of tree amplitudes to the D*^D^ and D*^D^ final states. If 5 = and equal amplitudes 
are assumed, one expects that ^ = 0, (7+ = C_ = and S+ = = — sin20i. 

The analysis described here is based on 140 fb^^ of data, corresponding to 152 x 10^ BB 
pairs, collected with the Belle detector at the KEKB asymmetric energy storage rings [3| • 
Two reconstruction techniques, full and partial, are used to increase the reconstruction 
efficiency. The event selection is similar to that in our previous publication j^, however 
some requirements are relaxed in order to increase the size of the sample used to extract 
the CP violation parameters. The full reconstruction method allows extraction of the signal 



decay with high purity, but, due to the small branching fractions of charmed meson decays 
into reconstructable final states, results in a low efficiency. In the partial reconstruction 
method, a D~ meson is fully reconstructed while only the slow pion (vt^q^) is required to be 
detected from the decay D*^ -^°^^ow- 

Neutral D mesons are reconstructed in five decay modes: K~ti'^ , K^ti^h^ti^ , K^tt^tt^, 
J'^s'vr+vr" and K^K~ j^. Charged D mesons are reconstructed via decays into K^ti^ti^ , 
Ksi^~ and K^K^-k^ . The selected combinations are fitted to a common vertex and a vertex 
quality requirement is applied to reduce combinatorial background. A ±15MeV/c^ interval 
around the nominal D mass is used to select D meson candidates for all modes except — >■ 
Ji^7r+7r°, for which ±24MeV/c^ is used (~ 3cr in each case). The selected charmed meson 
candidates are then subjected to mass-vertex constrained fits to improve their momentum 
and vertex resolution. We refer to a D candidate as having valid vertex reconstruction if it is 
formed by at least two tracks with hits in the silicon vertex detector. To suppress feed-down 
from the Cabibbo allowed decay D*^Di*^ due to K/n misidentification, we apply 

a D~ veto for the A'+tt^tt" and Ksti~ channels: if a pion candidate can also be 

identified as a kaon, and if, after reassignment of the kaon mass, the invariant mass of the 
combination is within ±15MeV/c^ of the nominal mass, the combination is rejected. 
This requirement suppresses the D*~^Di*^ background by a factor of 10 with signal 

efficiency of 98%. The D*~^ candidates are formed from D^tt^^^ combinations with invariant 
masses within ±2MeV/c^ of the nominal D*~^ mass. 

In the full reconstruction method we define B^ candidates as combinations of oppositely 
charged D*~^ and D~ candidates, where at least one of the D~ or the from the D*^ 
decay has valid vertex reconstruction. The signal is identified using the cm. system energy 
difference AE = E% — -Ebcam and the beam-energy constrained mass Mbc = V^^bcarn^^^^^' 
where E*^ (Pb) energy (momentum) of the B candidate in the cm. and -Ebeam is 

the cm. beam energy. B^ candidates are preselected by requiring \AE\ < 100 MeV and 
Mbc > 5.21GeV/c^. In case of multiple B^ candidates in this region a single candidate 
per event is selected based on the minimum sum of x^/DOF of the fits to intermediate 
charmed mesons. The scatter plot of AE vs. Mbc and the AE and Mbc projections are 
shown in Figs.^ In the Mbc projection 5° candidates are selected from the \ AE\ < 20 MeV 
region; the AE distribution is plotted for the region Mbc > 5.27GeV/c^. A fit to the My^ 
distribution with a Gaussian representing the signal contribution and an ARGUS function [2| 
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parameterizing the background finds a signal yield of 161 ± 16 events. A fit to the AE 
distribution is performed using a double Gaussian to parameterize the signal, while the 
background is described by a linear function. This fit yields 149 ±18 signal events. The cos^ 
distribution, where 6* is a decay angle in the D*^ rest frame relative to the boost direction, 
for the candidates from the \^E\ < 20MeV and Mbc > 5.27GeV/c^ region determined 
from Mbc fits is shown in Fig. ^ d) and is in good agreement with the Monte Carlo (MC) 
expectation. 
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FIG. 1: Kinematic distributions of D*^D^ candidates: a) scatter plot of AE vs. Afbc, b), 

c) AE and Mbc projections, d) cos 9 determined from Mbc fits in the data (points with error bars) 
and in the signal MC (histogram). The curves represent the fit described in the text. 

In the partial reconstruction analysis we define B candidates as combinations of D~ 
with valid vertex reconstruction and 7r^o„. As in our previous publication P], the angle a 
between the D~ and vt^q^ cm. momenta, and the D*~^ helicity angle 6, calculated using 
kinematic constraints, are used to identify the studied decay. We use the — > K~^7t~tt~ 
decay mode only. In addition, the D~ cm. momentum is required to lie in the inter- 
val 1.63GeV/c < < 1.97GeV/c, and the cm. momentum of vr^^^ is required to be 
smaller than 0.2GeV/c. Both momentum intervals correspond to the kinematic limits for 
the studied decay. In order to make the fully and partially reconstructed samples statis- 
tically independent, vr^^^ is rejected if, after being combined with any in the event, it 
forms a D*~^ candidate. The presence of a lepton (Aag) in the event is required to pro- 
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vide flavor tagging, suppress the continuum background to a negligible level, and also re- 
duce the combinatorial BB background. Charged tracks with cm. momenta in the range 
1.1 GeV/c < P/^^^ < 2.3GeV/c, which are identifled as muons or electrons are considered 
as leptons. A large fraction of the selected -D~£tag combinations originate from the decay of 
the same B meson: B D~£^X. This background is removed by a kinematic requirement 
that D~ and £tag do not originate from the same B: 



where E^^^^^ (-^o^ag) energy (momentum) of the D~ita.g combination. The effi- 

ciency of this requirement for the signal is estimated from MC simulation to be 87%, while 
the background is suppressed by a factor greater than 2. This requirement also removes 
leptons produced from the unreconstructed in the signal decay, an additional source of 
mistagging. We select partially reconstructed B^ candidates by requiring cos a < and 
|cos6'| < 1.05. In case of multiple candidates, the D^n^^^ combination with the best prob- 
ability of the D~ vertex fit or the largest |cos6'| is selected. The expected number of signal 
events in the partial reconstruction sample is calculated from the full reconstruction signal 
yield relying on the MC ratio of full and partial reconstruction efficiencies. We estimate 
A'partiai = 133 ± 13 and use this number to fix the signal fraction in later fits. 

We cross-check this result by estimating the signal fraction from the data. The distribu- 
tions of cos a for two regions of cos 6 are shown in Fig. |2l after imposing a tight requirement 
of ±8MeV/c2 (~ 2a) on the D~ mass. The first region 0.50 < |cos6'| < 1.05 (Fig. Ha)) is 
signal enriched due to the D*^ polarization; the second region |cos6'| < 0.50 (Fig. |21b)) is 
dominated by background. In a simultaneous fit to the two cos a distributions the signal 
shapes are fixed from the MC. The combinatorial background is parameterized by a second 
order polynomial function. The contributions from B^ D*^Di*^" and D*^D*~ are 

fixed from the MC simulation. The fit yields 137 ±39 signal events, in good agreement with 
the yield expected from the full reconstruction analysis. 

In the full reconstruction method, charged tracks that are not associated with the re- 
c„„.tn.cted ^ O-^D^ a.e used to ,de,.if. the flavo. of the s,g.a, B« Q. Eve.t. are 
divided into six subsamples of the parameter r, which is an event-by-event, MC-determined 
fiavor-tagging quality factor that ranges from r = for no fiavor discrimination to r = 1 for 
unambiguous fiavor assignment. The wrong-tag fraction and difference between B^ and B^ 



(-E'beam -^D£tag)^ ^b' ^1 
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FIG. 2: Distributions of cosa for: a) 0.50 < |cos^j < 1.05; b) |cos0| < 0.50. The fit functions are 
shown with sohd Hnes; the combinatorial backgrounds are presented by dashed hues. 

decays in each interval {wi and Awi, i = 1,6) are fixed using a data sample of self-tagged 
decay modes. In the partial reconstruction case, flavor tagging is provided by the high 
momentum lepton required in the event; the wrong tag fraction is determined from data as 
discussed below. 

The proper-time difference between the reconstructed and tagged B decay is calculated 
as At = {zo*D — Ztag)/P'yc, where zd*d and Ztag are the z coordinates of the two B decay 
vertices and f3'y = 0.425 is the Lorentz boost factor at KEKB. In the full reconstruction 
method we determine the B^ signal vertex by fitting the momentum vectors of D'^ and/or 

candidates with well-reconstructed vertices with the constraint of the interaction region 
profile. The tagging B vertex is found using well reconstructed charged tracks not assigned 
to the signal B^ and excluding tracks that form a Ks candidate. The signal resolution 
function parameters are obtained from the At fit to the 5° lifetime for the events from the 
signal region. The tagging B^ vertex resolution function is fixed from In the partial 
reconstruction method, the signal and tagging B^ vertices are reconstructed using the D~ 
candidate and £tag, respectively. In this case, both the resolution function parameters and 
the wrong tag fraction are extracted from the data using a sample of B^ D~l^vX decays 
tagged with a high momentum lepton. The D^i'^ combinations are required to originate 
from the same B decay based on the recoil mass against the D~i~^ system and its cm. 
momentum. The selected D~i'^ combinations are almost pure B^ D~l^v signal events 
with a small admixture of B^ [0^1x^)1^ the latter process is also considered as 

signal. A small contribution from combinatorial background under the D~ peak is estimated 
using D~ mass sidebands. The D~ , and 4ag vertices {zc-, Z£+ and -z^tag) are reconstructed 
using identical procedures to those used in the partial reconstruction method. The resolution 
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function is extracted from an unbinned maximum likelihood fit to the Ati = {z£,- — zi)/ (i^^c 
distribution. The wrong tag fraction is found from a fit to Attag = (^d- — ^tag)//57C to 
be w = 6.1 ± 0.9%. As a cross-check, the lifetime and mixing parameter Am^ are also 
measured from the fit to the Attag distribution to be 1.48 ± 0.04 ps and 0.52 ± 0.02 ps~^, 
respectively, consistent with lll |. 

In the full reconstruction method, the signal region is defined as \^E\ < 50MeV and 
Mbc > 5.27GeV/c^, and contains 360 events with 46% signal purity. In the partial 
reconstruction method, the signal region is chosen as \Mx+-k--k- — M£,-\ < 15MeV/c^, 
cos a < —0.9 and |cos^^| < 1.05. The total number of selected events is 2174 with 6% 
signal purity. 

We determine the CP violation parameters from an unbinned maximum likelihood fit to 
the At distribution. The signal probability density function is given by Eq. ^with effects due 
to mistagging taken into account. The resolution function Ro* d is formed by convolving four 
components: the detector resolutions for Z£,*£) and ^tag? the shift in the Ztag vertex position 
due to secondary tracks originating from charmed particle decays, and the smearing due to 
the kinematic approximation used to convert Az to At [l^. For each event we define the 
following likelihood value 

P^= [ \^^^^^VD'DiAt')RnM^U- At') 

J '-JD*D + Jhg 

+ r , Pbg(AtOi?bg(At. - At')] dAt', (4) 

JD*D -r Jbg -I 
where signal {fD*D) and background (/bg) fractions are calculated as functions of the follow- 
ing variables: AE and Mbc (full reconstruction); M/j- and cos a (partial reconstruction); 
cos^ (both cases). The signal distributions of the variables used for fo'D parametrization 
are determined from the MC simulation. The background parameters are obtained from the 
data. 

In the full reconstruction, the background At shape is fixed using the large Mbc-AE" 
region excluding the signal region. For the partial reconstruction the background con- 
tains a combinatorial component, for which the shape is obtained from D~ mass sidebands 
(30MeV/c^ < \Mk+-k-tt- — Md-\ < 60MeV/c^), and a component containing a real D~ , 
which may come from B decay. The shape of the latter includes a mixing term, and is 
obtained from a sideband (—0.8 < cos a < 0.0). 

Finally, the results for the CP violation parameters A, S±, and C± obtained from the 



TABLE I: Fit results. 



full rec. partial rec. combined 

A +0.03 ±0.09 +0.16 ±0.18 +0.07 + 0.08 

S- -1.17 + 0.48 -0.65 + 0.79 -0.96 + 0.43 

C_ +0.33 + 0.29 +0.11 + 0.45 +0.23 + 0.25 

S+ -0.25 + 0.52 -0.92 + 0.58 -0.55 + 0.39 

C+ -0.34 + 0.27 -0.39 + 0.38 -0.37 + 0.22 



individual fits to the statistically independent full reconstruction and partial reconstruction 
samples, as well as the result of the combined fit, are summarized in Table HJ We calculate 
the combined statistical significance oi CP violation to be 2.1a. The significance is defined 
as Y^'— 2Tn^Z^^7Anaxy) where Lmax is the likelihood returned by the combined fit and Lq 
is determined from a fit with the parameters A, S± and C± constrained to the values 
corresponding to no CP violation: ^ = 0, 5*+ = — S'_ and C+ = — C_. 

The At distributions for the subsamples having the best signal and tagging quality 
(/d*d > 0.1 and r > 0.5) after background subtraction are shown in Fig. El a) and b) 
for the full and partial reconstruction methods, respectively. 

The systematic error is dominated by the uncertainties in the signal fraction (±0.07 for S± 
and ±0.03 for C±), wrong tag fraction (±0.05 for S± and ±0.03 for resolution function 
parameterization (±0.05 for S± and ±0.02 for C±) and vertexing (±0.05 for S± and ±0.01 for 
C±). Other contributions come from the correlated backgrounds and signal box definition. 

n 

The result is consistent with We perform a number of cross-checks for our measurement. 
Using an ensemble of MC pseudo-experiments, we check both the linearity of the fitting 
procedure and the reliability of the statistical errors returned by the CP fit. A similar CP 
violation study is performed with self-tagged D*^Dj decay using both full and partial 

reconstruction techniques. The combined fit yields A = 0.00 ± 0.03, = +0.08 ± 0.12, 
C„ = -1.11±0.07, S+ = +0.00±0.12, andC+ = +1.12 + 0.07, consistent with the expected 
values A = 0, S± = and C+ = — C_ = 1. 

In summary, we have performed a search for the CP- violating asymmetry in the decay 
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FIG. 3: Background subtracted At distributions in the a) full and b) partial reconstruction meth- 
ods. The curves show the result of the fits. 

D*^D^ using two methods of reconstruction. From the combined fit to the data we 
have measured A = +0.07±0.08±0.04, 5_ = -0.96±0.43±0.12, C_ = +0.23±0.25±0.06, 
= -0.55 ± 0.39 ± 0.12 and C+ = -0.37 ± 0.22 ± 0.06. These are the most precise 
measurements of these parameters to date. The significance of nonzero CP violation in 
5° ^ D*^D^ is 2.7cr. 
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